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Abstract 

The effects of magnetic flux in statistical magnetisms, including Pauli paramagnetism, Landau 
diamagnetism, and De Hass-van Alphen oscillation, are discussed. It is shown that the diamag- 
netism could be much increased by the fractional magnetic flux, and the amplitude of the magnetic 
oscillation of De Hass-van Alphen can be amplified by the quantum effect of the flux. 
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I. INTRODUCTION 



Nano technologies has led to devices with enhanced functionalities under new operating 

principles such as quantum interferences that take place between the magnetic flux and 

electrons which in principle may effectively increase the operational velocity, and drastically 

decrease the power loss |l| . In this paper, we discuss the effect of magnetic flux in statistical 

magnetisms of electron gas. It is shown that the Landau diamagnetism (2), 3], and De Hass- 

van Alphen (dHvA) oscillation 0, Q could be much influenced by the fractional value of 

the flux that exists in the physical system has been confirmed by experiments in recent 
ft 

years |6|. Since the dHvA effect is an important technique in studying the energy bands 



of materials and the geometrical shapes of Fermi surfaces [5|, the present results may be 
applied to estimate the influence of topological defects or magnetic impurities on the Fermi 
surface. 

This paper is organized as follows. In Sec. II, the energy spectrum of an electron in 
a uniform magnetic field plus a magnetic flux is given. With which the partition function 
for the degenerate electron gas under the fields is calculated. In Sec. Ill, the statistical 
magnetisms of weak- and strong-degeneracy electron gas are presented. The influences of 
fractional magnetic flux in paramagnetism, diamagnetism, and dHvA oscillation are dis- 
cussed. Conclusions are summarized in Section IV. 



II. PARTITION FUNCTION FOR AN ELECTRON GAS IN THE PRESENCE OF 
A MAGNETIC FIELD AND A MAGNETIC FLUX 

The Hamilton operator H (x, p) of a spin-1/2 electron with effective mass /i and charge 
— e moving in an uniform magnetic filed B along the z-axis is given by 
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H (r, p) = — (pi + pj) + -^\{x 2 + y 2 ) + u L l z + — f z + fi B aB, (1) 

where the Larmor frequency ujl = eB/2fic, the Bohr magneton = \e\ h/2fic, o = ±1 are 
the eigenvalues of z-component spin operator a z , and the operator l z = xp y — yp x is the z- 
component of orbital angular momentum. To obtain the expression, we have used the axially 
symmetric vector potential A = (B x r)/2. From the Hamiltonian, we see that the charged 
particle continues to move freely in the ^-direction with a corresponding kinetic energy p 2 /2[i. 
The energy eigenstates \l/( )(r) can be chosen as ^W(r) = ^±(x, y) exp{ip z z/h}, where 
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^±{x,y) is the eigenstate of the charged particle in the x-y plane. In polar coordinates, the 
transverse state ^± can be decomposed as ^f±(p, (p) = R m k{p) exp{imip} in which R m k{p) is 
the 2D radial eigenstate and exp{imip} the angular part. The subscripts m — 0, ±1, ±2, • • •, 
and k = y/2pe/h are used to denote the different eigenstates in 2D plane. The energy 
eigenequation H^^(r) = i£\&(°)(r) then can be expressed as 
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For our purposes, it is convenient to write the wave function in a form 

oo 

tf(r)= E c mJ R m ,(p)e^e^^, c m = constant. 

m=— oo 

This is due to the fact that the system is linear. With \l/(r), Eq. <^fy becomes 
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So far, we have not included the influences of a magnetic flux in the electron. The nonlocal 
effect of_ajnagnetic flux is conveniently considered by a global nonintegrable phase factor 
The NPF represents the interaction of a charged particle with the magnetic 



(NPF) 



•i a umi 



field by the phase modulation as follows: 

ijjr(r) = M>( r ) exp 



i e 



he Jc 



A(r') • dr' , 



(5) 



where the wave function \I/(r) is the new state that has interacted with the magnetic flux 
defined by the vector potential A(r), and the subscript C indicates the value of the phase 
depends on the choice of different paths. For an infinitely thin tube of a finite magnetic flux 
along the z-direction, the vector potential can be expressed as ([101]. Ch. 15) 
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Here e x ,e y stand for the unit vector along the x, y axis respectively. Introducing the az- 
imuthal angle ip(r) = tan _1 (?//x) around the magnetic tube, the components of the vector 
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potential can be in terms of Ai = 2gdiip(r). The associated magnetic field lines are thus 
confined to an infinitely thin tube along the z-axis 



B 3 = 2ge 3ij d i d j <p(r) = 4ng5(r ± ), (7) 

where r_i_ represents the transverse vector rj_ = (x,y). Since the magnetic flux through 
the tube is defined by the integral $ = J dxdyB 3 , the coupling constant g is related to the 
magnetic flux by g = $/47r. By using the expression of A{ = 2gdiip, the angular difference 
between an initial point r' and the final point r in the exponent of the NPF is given by 

ft . r* -yx + xy f*rxdr 
(p-tp = rfrv?(r) = / dr = / (8) 

Jt' Jt' x z + y z Jx> r z 

where ip = dip /dr. Without losing the generality, we choose <p' = 0. Given two paths C\ 
and C*2 connecting r' and r, the integral differs by an integer multiple of 27r. The winding 
number is given by the contour integral over the closed difference path C: 

1 r rxdr 

The interaction of electron with the magnetic flux is therefore purely nonlocal and topologi- 
cal. Its action in the NPF takes the form A mSug = —a 2irn, where a = —2eg = — $/$ 
is a dimensionless number with the customarily minus sign. The NPF now becomes 
exp {— iao(2nn + ip)}. The wave function *& n (p, <p, z) for a specific winding number n can 
be obtained by converting the summation over m in (j4j) into an integral over a and another 
summation over n by the Poisson's summation formula ([10], Ch.2) 

oo roo oo 

£ f(m) = I da e 2 ^f(a). (10) 

m=— oo 00 n=— oo 

Equation (jlj) can then be cast into 

da ^2 °a \ \E — ahui — 

^ n=-oo l\ ^/i / 
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x R ak (p)e t{a - ao) ^ +2n7T) e iPzZ/h = (11) 



Obviously, the number n in the right-hand side is precisely the winding number by which we 
want to classify the wave functions. Employing the Poisson's formula J2n exp{zfc(<^+2n7r)} = 



J2m=-oo — m)exp{imip}, the summation over all indices n forces a = a modulo an 
arbitrary integer number, and yields 



+ 
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2/j, \dp 2 pdp p 2 

x ^.t^e^e^ = 0. (12) 

One see that the effect of the magnetic flux in the wave function is to replace the integer 
quantum number to with a real one (to + ao) that depends on the magnitude of magnetic 
flux. With the orthonormal integration J 27r e^™ - ™'^ dtp = 2ir5 mm ', the corresponding radial 
wave equation is found to be 

| - (to + a )hu L - |^ - VBcrB^J 
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The energy spectrum is given by 

v 2 

E = [2n p + \m + a | + (to + a ) + 1] acj L + ^ + // B o"B (14) 

with n p = 0, 1, 2, • • •. When we restrict the problem to the 2D x-y plane, the spectrum 
reduces to that of Landau levels if the flux is quantized at integer value. Along the direction 
of the magnetic field, the electron is free motion. With the help of the spectrum, the partition 
function of the electron gas under an uniform magnetic field and a magnetic flux can be 
calculated. By defining the cyclotron frequency u c = 2u>l and expressing (m + a ) = (rh + Q 
with the possible values to = 0, ±1, ±2, • • • and < £ < 1, the spectrum is expressed as 
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E no = (n + C + -)hu c + g + p B aB, (15) 

where no = {n p + to) with the range no = 0, 1, 2, • • •. It is (, the topological nonlocal effect 
of a fractional magnetic flux, leads to the amplification of the Landau diamagnetism and 
the dHvA effect. 

As usual, we start to discuss the magnetisms of the electron gas by the grand partition 
function 

oo 

\nZ= ^ln{l + exp[/3(r-K )]}, (16) 

no=0 



where 1/(3 = kT is the mean thermal energy, r is the chemical potential, and E no is the en- 
ergy spectrum of a single electron. To calculate the partition function In Z, it is conveniently 
to use the Mellin transformation [ill . [l^ | 



1 rc+ioo 

f(x) = — / Futyx'dt, (17) 
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and its transformation pair 

fix) F M (t) 



ln|l + a;| <-> -csc(Trt), -1 < Ret < 0. (18) 
The grand partition function can be expressed as 

1 00 rc+ioo 1 , s ,7+ 

XnZ = k^ ■ id^ e t' 0<C<L (19) 

In which X^=o exp{ — (3E no b) = Z\{(3t) is the canonical partition function for single particle 
in Boltzmann statistics. The grand partition function turns into the representation 

1 rc+ioo I rlf 

\nZ = - ——e^Zxipt)^ °< c<1 - (2°) 

2% Jc-ioo sm(7rt) t 

With this, the partition function of a system can be obtained if we can perform the integra- 
tion. Since the degeneracy of each Landau level g = (fiV 2 ^ 3 uj c )/h, the canonical partition 
function Z\{(5) of single particle is found to be 



no=0 it=±1 



f 1 0P 2 1 

x exp | -/3(no + C + -)^ c - ^ - /fys<7B j 

_ 2K / Q cosh/ 2do 

" A3 sinh/ 6 • [Zi) 
Here the volume of system V = L x L y L z , the thermal wave length A = (h 2 /2^^ikT) l l 2 , and 
l Q = (H/3u) c )/2 = (3^bB as well as l s = (3[XbB indicates that they come from the orbital 
and spin motion respectively. One note that in the limit B — > 0, we have Zx(J3) — > (2VQ/A 3 
which is a well-known result. The grand partition function now becomes 

W /• c+i0 ° cosh^tW^- 2 ^' , , , 

lnZ = — 1/ — . v 7 \ . - - . dt, 0<c<l. 22 
%X 3 Jc-ioo t 3 / 2 sin 7rt sinh (l a t) ' v ; 
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The integrand have single poles at t — ±n from sin(7rt) and t = ±(inp)/l from sinh(/ t). 
Furthermore, a branch point locates at the origin due to t 3 ! 2 . For the case of weak degeneracy, 
i.e. (@t — 2ql Q ) < 0, the integral can be performed by chosen a large closed contour along 
the complex t-plane of the right hand side. It yields 

cosh(n/ s )e (/3T - 2do)n 



71=1 



n 3 / 2 sinh(n/ Q ) 
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when £> ^ 0. Since electron just has a tiny mass, the electron gas is in general strong 
degeneracy. So we will stop here for the case of weak degeneracy. For the case of strong 
degeneracy, i.e. (f3r — 2ql Q ) > 0, one adopt a closed contour along the left complex plane as 
shown in Fig. 1. The residues arising from poles t = ±(iirp)/l , p — 1, 2, • • •, are given by 
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The contour integral can then be written as 
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where 7 indicates that the range of the integration is along the branch cut from C to D. 
Since (f3r — 2ql ) > 0, the integrals J BC and J DA vanish as the radius r — > 00. The grand 
partition function becomes 
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Again, due to the fact that {(3t — 2ql Q ) > 0, the dominant contribution of the contour integral 
J is around the neighborhood of \t\ =0 such that the integral 

l I cosh(/ s t)e( /3r - 2d °)* 
~2i 7 7 t 3 / 2 sin(Trt) sinh(U) 
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Here y = (/3r — 2<^ G ), and we have used the representation of Hankel's contour integral for 
Gamma function 

~ I t- z e l dt. (29) 
h 



1 



r(z) 2m 

The grand partition function for the gas of spin-1/2 electrons in magnetic field and flux is 
finally given by 
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The last term is an oscillatory term which becomes significance when the system is presented 
in a strong magnetic field (l s — l Q > 1). In the following section, we shall see that the 
paramagnetism and diamagnetism come from the third term of the first middle bracket that 
turns into dominant when the system is presented in a weak magnetic field (l s , l « 1). 
The first two terms are the result of a electron gas with strong degeneracy. 



III. EFFECTS OF A MAGNETIC FLUX IN PARAMAGNETISM, DIAMAG- 
NETISM, AND DHVA OSCILLATION 



We first discuss the conditions of strong degeneracy (y 1) and weak magnetic field 
(l s ,l << !)• In this case, the oscillatory term can be neglected, and 
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With this approximation, the density of the degenerate electron gas is related to the Fermi 
energy as follows: 
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where e/ = (3n 2 nh 3 ) 2//3 /2[i is the Fermi energy of electrons at zero temperature T = 0. Eq. 
(1321) gives the value of (3t in terms of fief 

r 2 



/3t « (Pe f + 2sl ) 



^-{(3e f + 2,l )- 2 



(33) 



With this connection, one can connect the flux effects with the magnetisms of electron gas. 
As usual, the magnetization M and the magnetic susceptibility % are expressed as follows: 

/i d 



M= V (3 ~ 1 ^B lnZ > X = /io f (3VB8B 



lnZ, 
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where [Xq is the permeability of the vacuum. With the conditions y ^> 1 and l s , l a << 1, the 
dominant contribution of fl30l) comes from the first middle bracket. It follows that 
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where the main contributions come from the terms including y 1 / 2 and y 3 ^ 2 ■ With the help of 
( 1331) . the term with variable y 1 ^ 2 can be identified with the well-known Pauli paramagnetism 
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and the Landau diamagnetism 
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The contribution of the term including y 3 ! 2 gives 

2n/i /i B 



Xab 



LL B B 



IT 

1 



D' 



(38) 



which is an effect of diamagnetism. The subscript AB indicates the contribution comes 
from the magnetic flux of Aharonov-Bohm. Under the general condition, ^bB -C kT C e/, 
Xab may become more important than the Landau diamagnetism xl when the value of the 
nonquantized flux q — > 1. However, it will vanish if ^ — > 0. We see the nonlocal effect 
of flux in the statistical magnetism of an electron gas is significant, and depends on the 
nonquantized value of the flux. On the other hand, the paramagnetism does not affected by 
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the flux which is reasonable since the flux does affect the space degree of freedom, not the 
internal degree of freedom. 

In the condition of an strong magnetic field, HbB ~ kT <ej, the oscillatory part 
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becomes important, where l = l s is considered. The dominant contribution to the suscep 
tibility is the cosine function. So the susceptibility reads 
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where we have used (!33|) to approximate the {3t and omits the small quantity 
—7rp(kT) 2 /(12fiBBef) in the sine function. The second term in the former factor is the 
effect of the flux which linearly increases the amplitude of the dHvA oscillation. In fact, 
we can show that when <^ — *> 1 the flux effect is more important than contributions from 
Z 3 / 2 and [sinh (tt 2 p/1 )] 1 in fl39l . Indeed, the differentiation d/dB with respect to l z J 2 and 
[sinh (ttVOP in (JH give the contributions 
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Both are smaller than the second term in (14D1 when ^ — > 1. Before finalizing the paper, let 
us compare the magnitude of magnetization between the oscillatory part and nonoscillatory 
part. The magnetization of the later is given by 
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The oscillatory part approximates to 
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The ratio between these two parts has the order of magnitude 
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When B = 10 Tesla, the quantity /Ub-B ~ 5.8 x 10 eV. The Fermi energy tf of the electron 
gas in metals is of the order of a few electron volts. The ratio M osc /M is a significant 
quantity when the value of fractional magnetic flux ( — > 0. In this case the oscillatory part 
dominates the magnetization. On the other side <j — > 1, nonoscillatory part M would be 
dominant the magnetization. 



IV. CONCLUSIONS 



In this paper, the effects of a magnetic flux in statistical magnetisms, including Pauli 
paramagnetism, Landau diamagnetism, and the dHvA oscillation, of a degenerate electron 
gas in a magnetic field of arbitrary strength are studied. It is shown that the diamagnetism 
can be increased by the fractional magnetic flux, and the amplitude of the magnetic oscilla- 
tion of dHvA can be amplified by the quantum effect of the flux. Since the effect of vortex 
in a degenerate Fermi gas can be realized in the experiment nowadays [l^] and is important 
in exploring the critical states of matter, the result may be useful in studying the statistical 
properties of a degenerate Fermi gas involving the vortex. 
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FIG. 1: The contour integral on the complex i-plane for the strong degeneracy. The poles at t = n, 
t = —n, and t = ±(mp)/l , p = 1, 2, • • • are shown with different notations. 
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